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Abstract. We develop a method to constrain non-isotropic features of Cosmic 
Microwave Background (CMB) polarization, of a type expected to arise in some 
models describing quantum gravity effects on light propagation. We describe the 
expected signatures of this kind of anomalous light propagation on CMB photons, 
showing that it will produce a non-isotropic birefringence effect, i.e. a rotation of 
the CMB polarization direction whose observed amount depends in a peculiar way on 
the observation direction. We also show that the sensitivity levels expected for CMB 
polarization studies by the Planck satellite are sufficient for testing these effects if, 
as assumed in the quantum-gravity literature, their magnitude is set by the minute 
Planck length. 
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1. Introduction 

One of the main open problems in modern physics is finding a description for phenomena 
taking place in the so-called "quantum-gravity realm" , at scales where both gravitational 
and quantum effects are non- negligible [HEj- Several alternative candidate theories are 
being developed but discriminating among them experimentally is extremely difficult 
as a result of the smallness of the characteristic quantum-gravity length scale, the 
Planck length Lp ~ 10~'^^m (which of course corresponds to an energy scale, the Planck 
scale Ep ~ lO^^eV, which is huge by particle-physics standards). Nonetheless it is 
now well established that astrophysical observations can provide significant information 
about physics at the Planck scale [3], and recently it has been pointed out that also 
cosmological observations can be helpful. The opportunities on the cosmology side are 
becoming more mature as we start to have high precision data, and in particular it was 
recently shown that current data on Cosmic Microwave Background (CMB) radiation 
have already enough sensitivity to be used to establish meaningful constraints on the 
effects predicted by some Planck-scale theories [H O Ej . 

This is possible because, even if CMB radiation is characterized by quite low 
energies, it has been propagating for very long times, so that it is potentially subject 
to a large accumulation of new physics effects affecting photons propagation. In [1] 
some of us were able to show that measurable changes can be induced on CMB 
polarization power spectra due to anomalous light propagation, of the kind predicted 
by some descriptions of modified electrodynamics [3, [H] which are much in use in the 
literature as effective description of the implications for electrodynamics of the Planck- 
scale (possibly quantum) structure of spacetime. These anomalous properties of light 
lead to a birefringent behavior of CMB radiation, producing a rotation of its polarization 
direction, which is detectable through an analysis of the cross-correlation power spectra 

p mini [121 [13]. 

In [1] anomalous light propagation was encoded in the deformed-electrodynamics 
model introduced by Myers and Pospelov [7], a dimension- five effective field theory 
coupling the electromagnetic field with an external fixed four-vector, which explicitly 
breaks Lorentz symmetries. This model was also tested in several astrophysical contexts, 
leading to very stringent constraints [HI [iSl [HI [lEl [IZj • In [H] it was observed that these 
bounds exploit significantly the spatial isotropy regained by the ad hoc choice of having a 
purely timelike symmetry- violating four- vector, and actually this choice is only available 
for a restricted class of frames of reference, since the four- vector will of course still acquire 
a spatial component in other boosted frames. This assumption can be limiting in two 
ways. Constraints on one single component of a four vector derived in different reference 
frames (like the rest frame of two different astrophysical sources) cannot be compared 
without any information on the other components of the four-vector. Moreover in [18] 
it was shown that the limits placed assuming the four-vector to be purely timelike do 
not give reliable information on the most general case in which the four-vector has all 
the components different from zero, so that also space isotropy is violated. 
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The success of the attempts to constrain isotropic anomalous hght propagation 
through CMB observations provide motivations to investigate more widely how Planck- 
scale modifications of electrodynamics can affect propagation of light in a way that 
would be detectable through cosmological observations. 

In this paper we describe the expected behavior of CMB photons in presence of 
the non-isotropic Lorentz-symmetry violations described by the most general version 
of the model proposed in Ref. [7]. We show that anomalies can affect CMB radiation 
propagation, in a way that can be described as non-isotropic birefringence, i.e. a rotation 
of the polarization direction whose amount is characterized by a peculiar dependence 
on the observation direction. 

The most relevant characteristic of the model described in this paper is that 
anomalous light behavior depends on its propagation direction with respect to a 
preferred direction codified within the symmetry-breaking vector tIq. To make more 
explicit the difference with the previously studied case [TJ H] in which space isotropy 
is preserved (n^ = (rio, 0, 0, 0)), we will specialize at the end of section 2 to the 
complementary case in which the symmetry breaking vector has only the space 
components different from zero. 

We also show that in this case the peculiar space-dependence of the birefringence 
effect is such that the rotation can not be detected exploiting the standard tool of 
correlation power spectra, which instead was shown to be valuable for constraining the 
isotropic rotation effect studied in 

So we need to develop a specific method to detect such anisotropic effects, and this 
is described in section 3. In section 4 we forecast on the sensitivity level that will be 
reachable with the Planck satellite observations, showing that it will allow to constrain 
the model even beyond the Planck scale level. 

2. CMB photons propagation in presence of non-isotropic Lorentz 
symmetry violations 

In [18] it is described the behavior of photons in presence of non-isotropic violations of 
Lorentz invariance, resulting from the generalization of the Myers-Pospelov model, with 
nonzero spatial components of the four-vector. The Lagrangian density of the model is: 



where Ua is the symmetry-breaking four-vector and the coupling constant between the 
electromagnetic field and the vector is given by the inverse of the Planck energy, so 
setting the appearance of new physics at that scale. 

The relevant equations for the dispersion relation and the eigenstates of the field 
propagation resulting from the above Lagrangian are [18] : 
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where the field eigenstates are written in the basis 
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such that the first component of the field is the longitudinal one and the other two 
components describe the field in the plane transverse to the propagation direction. 

It is then clear that the model leads to the emergence of several anomalous 
effects, whose intensity is always depending on the relative orientation between the 
symmetry-breaking four vector and the field wave- vector. Eq.(l3]) shows that the 
field eigenstates are elliptical, and moreover are not transverse, having a Planck-scale- 
suppressed longitudinal component. Eq.([2]) describes a modification of the dispersion 
relation, of opposite sign for the orthogonal polarization states of the field £± (this is a 
birefringent behavior), and modulated by the angle between the spatial part of and 
the field propagation direction p. 

All these anomalous properties are of course strongly suppressed by the Planck 
scale, but in the following we will show that some of them can be indeed amplified by 
very long propagation times, compensating for the Planck-scale suppression, so that 
they can produce observably large effects. 



2.1. Effective behavior of CMB radiation 

CMB photons have a partial linear polarization due to their Thomson scattering with 
electrons at the last scattering surface in presence of a quadrupolar anisotropy. After 
this last scattering they propagate (almost) freely toward us. 

To study photon propagation described by Eqs. ([2]) and ([3]) we would like to write 
a plane wave with the known polarization properties of CMB as a linear combination 
of the eigenstates £± of Eq. (|3]). But a transverse linearly polarized field cannot be 
written as a linear combination of these two eigenstates, meaning that it does not belong 
to the solutions of the modified Maxwell equations following from the non-isotropic 
generalization of the Myers-Pospelov model. To overcome this difficulty, taking into 
account the fact that what we actually observe within the sensitivities presently available 



X We use the notation of Jones three-dimensional vectors 



the polarization state of the field 
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is that CMB is transverse and linearly polarized, we will assume that the original field 
is a solution of the modified Maxwell equations, which can be written as an expansion 
in powers of and at zeroth order can be described as a transverse linearly polarized 
state. So the field will be of the form 

E = + ^i(i) (6) 

where E^^^ is transverse and linearly polarized and E^^^ must be such that ^ is a solution 
of the modified Maxwell equations (to the first order in i.e. the polarization state 

S of E can be written as a linear combination of the eigenstates (IHj): 

1 



E, 



-£^^^ = AS. + B£_ 



(7) 



Of course also the complex coefficients A and B can be expanded as a series in powers 
of ^. To the first order: A = + ^A^^\B = + ^B'-'^I 

To enforce Eq. ([7]) and constrain the form of 8^^\ 8^^\ A and -B, it is convenient to 
rewrite also the eigenstates in Eq. (jS]) pointing out the expansion in powers of -J- (they 
were already written only up to the first order in -^): 
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We set t = at the last scattering surface. At this time, the zeroth order of the CMB 
field in Eq. (jS]) is a linearly polarized transverse field: 



E(°)(f,t = 0) = 3fJ{ 



/ \ 



with ai and 02 real coefficients (we are still using the basis (jlj)). Asking equation Eq. 
(I7j) to be satisfied to the zeroth order in -g-, we find 



A(o) 

5(0) 



0.2 —^Qj 

a2+iai 
V2 



(9) 



which is indeed the relation between the components of a linearly polarized field written 
in a linear basis and the components of the same field written in a circularly polarized 
basis. 

For what concerns the first order in we will see now that it is not necessary to 

Hip 

calculate A^^^ and B^^\ nor to give the explicit form of £^^\ After propagating for a 
time t the CMB field takes the form: 



E{r, t) = 3ft{ 
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(10) 
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where u± are defined in Eq. 
witli uq 



We expand also u± in power of 



UJ± = (jJq ± 5uj, 



\p\ and 5ijJ = "^IpI^ (^o + ^ ) , so tliat tlie above expression for tlie field 



can be written up to the first order in ^ as: 
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So up to the first order the field will be a combination of a term with the same 
polarization state £ as the field at time t = 0, propagating with classical frequency 
Uq, and a term that describes the mixing of the zeroth-order components of the field. 
The polarization direction of the zeroth-order field is rotated during propagation, as it 
happens in presence of birefringence. 

So we find that the first-order correction to the field, E^^\ is not significant within 
the orders of magnitude of our interest, and so we can consider the physical system as 
a linearly polarized field whose polarization direction rotates during propagation with 
angular velocity 6u. In fact all the corrections to the zeroth-order field are of order as 
Su, but the magnitude of the terms with 6uj can be amplified by the propagation time t, 
and this is indeed the case for CMB photons. So all the corrections to the zeroth-order 
field are much less important than the corrections due to rotation of the zeroth-order 
linear polarization, and we expect that, at least within the sensitivities presently reached 
by the instruments, we can at most see only the rotation of the zeroth-order linearly 
polarized transverse field: 
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where we have defined Tq = p ■ f — uqI. The rotation angle a 
rotations from the 2 direction to the 1 direction. 

Another interesting feature emerging from this approximation is that in particular 
the longitudinal component of the field (which is univocally determined by imposing 
that Eq. ([71)) is constant with time, and is not amplified by the long propagation time. 



5ujt is positive for 



2.2. Anisotropic rotation of CMB Stokes parameters 

The relevant quantities that we can measure about CMB polarization are the Stokes 
parameters Q and U related to photons linear polarization properties. Since they depend 
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on the choice of reference frame, in order to describe CMB Stokes parameters all over 
the sky, we have to set a conventional all-sky reference frame. This is of particular 
importance for us, since also the linear polarization rotation angle induced by (non- 
isotropic) birefringence can take opposite conventional signs depending on the reference 
frame choice. 

In the previous subsection we have conventionally set the rotation angle a{t, Uq, n ■ 
p, \p\) to be positive if the rotation direction goes from the axis 2 = to the 

axis 1 = — r==4^= for a field coming toward us. 

It is clear that, since we do not know the direction of n, we cannot trivially change 
reference frame to the one conventional for CMB to the one, n-dependent, that we have 
used up to now to describe non-isotropic photons propagation. But actually what really 
matters for the determination of the sign of a is the relative handedness of the two 
different reference frames, and the conventional reference frame of Eq. that we have 
used to describe non-isotropic birefringence has an handedness which does not depend 
on the direction of n. 

The reference frame in which conventionally the CMB Stokes parameters are defined 
at each point on the sky is a right handed coordinate system with z pointing outward 
toward the sky and x,y directions defined as in appendix C.4 of [T9] . 

In each point of the sky a rotation which in the reference frame ([1]) {1,2} = 
{ , =, j used here to describe the transverse component of the electric 

field is from the 2 direction to the 1 direction (so that a is positive), will correspond in 
the CMB reference frame to a rotation from the y direction to the x one. 

In the reference frame on the sky such that the observation direction is defined by 
the polar angles {6,(f)}, the direction of the vector n is given by the angles {6n,4>n}- 
Then in the associated cartesian reference frame n can be expressed as 

h = (cos(0„) sin(6'„), sin(</)„) sin(6'„), cos(6'„)). (14) 

An analogous formula will hold for the photon propagation direction p, which is opposite 
to the observation direction {9,(f)}, so 

p = {— cos(0) sin(6'), — sin(0) sin(6'), — cos(6')) (15) 

Then the dependence of the rotation angle a = 5ujt = ^'^o + t on the 

observation direction takes the forn|§|: 

a{9, 0) = ttIpI^ ('^o — 1^1 (sin 6* sin 9n cos(0 — 0„,) + cos6' cos6'„))'^ t. (16) 
Ep 

So at each point in the sky {6, 0}, non-isotropic birefringence induced by the Lorentz 
violating Lagrangian of Eq. ([1]) that we are studying produces a mixing between the 

§ Note that here we are not taking into account the effect of photons energy redshift due to propagation 
in an expanding universe. To account for redshift the expression for the rotation angle has to be changed 
in a way similar to the ones described in [4] 
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Stokes parameters, whose amount, for each given n, depends on the observation direction 
and at the first order in the rotation angle a is given by: 

g'(^,0) = g(^,0) + 2a(^,0)f/(^,0) 

0) = u{e, 0) - 2a(^, 0)g(^, 0) (ir) 

Note that a positive rotation angle a here corresponds to a rotation from the y to the 
X direction. 

2.3. Consequences on the polarization harmonic coefficients 

In the case of isotropic birefringence the uniform rotation of Stokes parameters can be 
detected looking at the cross-correlation and auto correlation power spectra, since it 
produces a peculiar mixing between them |31 |9l HOl ITT] - 

In this subsection we will work out the modification to the full-sky power spectra 
induced by the direction-dependent rotation described above in this section, following 
the scheme presented in [201 EB [22], but specializing the calculations to the case of 
our interest. This will show that the kind of birefringence with anisotropic behavior 
described by Eq. ( fT6l) cannot be detected using the standard tool of full-sky power 
spectra, so that it will be necessary to develop an alternative method of analysis for 
polarization data, which is done in the next section. 

The rotation angle a is a scalar function defined on the sphere, so it can be expanded 



into spherical harmonics: 

a(^,0) = 5^aLM>LAf(^,0) (18) 
LM 

In the case of our interest, in which a takes the form ( !T6|) . we need to sum only up to 
L = 3, with the expansion coefficients given by (we have defined A = ^-^t)'- 

aoo =2Ano{nl + \n\'^)^/^ (19) 

aio = —A \n\ (5^0 + |n|^)\/37rcos6'.„ (20) 

ai±i = ±\A\n\ (5^2 + Inne^^^-Veyrsin^n (21) 

"20 =^no|npy^(l + 3cos2^„) (22) 

a2T2 = Ano \n\^e^''^"^sme^^ (23) 

a2=Fi = ±^no|npe=^*'^"'y^sin6'„cos6'n (24) 
1 [tt 

a-so = -—A\n\^J - {3cos9n + 5cos39n) (25) 

a3T3 = T^lnl^e^^^^-^-^^sin^^^^ (26) 

^3^2 = -A \n\^e^''"^- J — sin 9n^ cos On (27) 



9 




ag^i = ^ \nfe^'^- W y (3 + 5 cos 2^„) sin On (28) 

We write the modification of multipoles of electric and magnetic modes of 
polarization taking only the first order in a and considering the possibility of having 
primordial non-zero magnetic modes. 

Following the scheme presented in [201 El] we are able to evaluate the modification 
of the power spectra due to the polarization rotation: 

{Cf^y = - 2£_^i X] "^^^ (^^^ " ^f^) ^u^emem (29) 

m LM 

2£+L=even 

icn' = cr (30) 

{Cj^y = ^ X ^ oilmCJ^ H^(^^^^^^ (31) 

m LM 

2l+L=even 

{cry = cr (32) 
{cry = cr (33) 

where the coefficients if^ and the Wigner 3j-symbols CiJfi^ are defined in [201 EI] • Notice 
that these expressions reduce to the ones for the isotropic rotation case if aoo = «o and 
o^LM = for all other L, M. In this case the mixing between the correlation power 
spectra are given by 

{cry= -2a,{cr-cr) (34) 

{cry = cr (35) 

{cry = 2aocr (36) 

{cry = cr (37) 

{cry = cr {^ 

Note that in this case a positive corresponds to a rotation from the x axis to 
the y axis. 

In the general case described by Eqs. fl29|) - fl33|) . all the correction terms for the 
spectra involve sums over 2i + L = even, so take contributions only from aiM with even 
L. Moreover from the selection rules of the Wigner 3j-symbols, is different from 

zero only if M = 0. So the only multipoles of a that contribute to the rotation of the 
power spectra are aLo, with L even. 

In our specific case this means that only the multipole coefficient aiM with 
{L,M} = {0,0}, {2,0} contribute to modify the full-sky power spectra. In particular, 
if the Lorentz-violating Lagrangian of Eq. ([1]) contains a purely space-like vector, 
Ua = {0,n), which is the case that is here of primary interest, then the rotation angle 
( IT6l) reduces to 

a{9, (j)) = amax (sin 9 sin 9n cos(</) — </>„) + cos 9 cos 6'„)^ (39) 
amax = - -^\p\'^\n\H 
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and from Eq. (1251) . setting uq = 0, we see that it needs only L = 1,3 multipoles for the 
spherical harmonics expansion, so leaving no traces into the power spectra. An intuitive 
way to understand why this happens is to realize that, for Ua purely space-like, the 
rotation angle a takes opposite signs in opposite directions in the sky. Essentially, the 
modifications to the power spectra coming from one hemisphere of the sky compensate 
the modifications, of opposite sign, coming from the other hemisphere, and full-sky 
spectra are not affected by the polarization direction rotation. We concentrate on 
Lorentz violations induced by a purely-spacelike vector coupled with the electromagnetic 
field, since the standard analysis of full sky power spectra, that turned out to be useful in 
constraining the isotropic version of the model |1] , cannot be exploited. In the following 
section we introduce a new method of analysis, which indeed exploits the peculiarities 
of the model here of interest, and we apply it on simulated data compatible with the 
ones expected from the Planck satellite observation. 



3. Method of analysis 

The free parameters of the model that we want to constrain are the module and direction 
of the symmetry-breaking vector n. Its module defines amax, the amount of rotation 
displayed by a photon coming from the direction of maximum effect 0„} (in the 
antipodal direction the rotation is also maximum but has opposite sign). In order to 
avoid the cancellation of the effect that would arise when working with (nearly) full sky 
CMB maps, we may work separately on the two hemispheres of the sky. However we 
don't know the direction n that would provide the natural division of the sky. A crucial 
observation to overcome this problem is that if we consider a small circular region of 
the sky centered at then the power spectra estimated on the region are rotated 

according to Eqs. (1M1) -(!381) by an angl^: 

ao = Ax Umax, (40) 

where A is given by the average of the function (sin ^ sin 6*^ cos (0 — 0„) + cos^cos6'„)^ 
over the region. For instance, if the angular radius of the region is 20° we obtain 
A = 0.913, while 10° sets A = 0.978. By taking this correction factor into account, the 
laws of transformation of the spectra at first order in the rotation angle are {C are the 
unrotated spectra): 



Cf^ = 2ao{Cf''-Cf''). (41) 

All the other spectra {Cf^, Cf^, Cf^ and Cj^) remain unchanged at linear order. Note 
that for tto = we expect Cj^ and Cf^ to be zero. We exploit this feature to constrain 

II Here and in the following we indicate with ao the effective rotation under the assumption of constant 
rotation, with OLmax the maximum amount of rotation due to the anisotropic birefringence effect, and 
with the same symbols with a ( ~ ) the corresponding estimated quantities 
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direction and amplitude of n. To this purpose we generate a set of 1000 masks to 
select, out of a CMB map, disks of angular radius 20°l3ll ^'^^ centers {6'c, 0c} distributed 
randomly in one half of the sky {9^ G [0°,90°], 0c G [0°,360°]). In correspondence of 
each of these disks there is another one centered in the opposite direction. In this way 
we cover the entire available sky allowing for a certain degree of superposition among 
the disks. 

We want to test the hypothesis that a disk is centered in the direction of maximum 
rotation. For this purpose we estimate the polarized power spectra out of each of them. 
We employ an estimator based on the pseudo spectra formalism [231 123], properly taking 
into account the effect of instrumental noise and incomplete sky coverage. According to 
Eq. f l4T]) . for each disk in one hemisphere we define the quantities: 

Dj^{ao)=Cr -2a,Cr (42) 

using the measured power spectra {Cf^), under the assumption that the measured 
EE, TE and BB spectra are unchanged at first order in the rotation angle (see Eqs. 
( l4T|) ). Since in the opposite hemisphere the effect is expected to be identical except for 
a sign flip in ao, the estimators Dj^ and Df^ for the antipodal disks can be easily 
built just changing the sign of ao in Eq. ( 142|) . In order to determine the direction of 
maximum rotation we compute, for each region i [i = 1, 1000) and its antipodal one, 
the ao = aoi which minimizes the x^(q^o) given by: 

X^(ao) = Yl ^rs- + Y: D^^tsM'^r + (43) 

. r)TBy-l f=)TB , f=)EBy-l f=)EB 

where TixY/e {^XY/e) are the covariance matrices of the Df^ {Df^) estimators. We 
end up with an estimate of ao, i.e. aoi, for each region and we project these values 
onto a sky map, taking the average of the aoj values associated to overlapping regions. 
Our estimate of amaxi i-e. otmax-, is the maximum value of the ao map rescaled by the 
correction factor A as given in Eq. ( HOl) . The uncertainties on cxmax can be assessed 
via Monte Carlo (MC) simulations, as it will be explained in detail in the next section. 
The final step of the analysis is to derive an estimate of the direction of the symmetry 
breaking vector, and the confidence intervals associated to it. By testing the 

method on simulations, it turns out that we cannot evaluate simultaneously Umax and 
{^n,0n} by minimizing a joint ^^(tto,^',^), due to the degeneracies between the three 
parameters. For this reason we take the direction of the maximum value of the ao map 
as the best-fit of {6'„, 0.„} and we derive the uncertainty on this estimate by slicing the 
function ao(6', 0) in correspondence of the 1 and 2(T errors on ao- 

^ Note that we choose this value as a trade off between the goodness of the approximation in Eq. PO)) 
and the sample variance increase due to power spectrum estimation on very small regions of the sky. 
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4. Forecasts for Planck 

In this section we derive forecasts to constrain the anisotropic birefringence model 
set forth in Section [2] with Planck [25], by applying the method outlined above on 
simulated data. In particular, we draw (T,Q,U) CMB simulated maps as Gaussian 
realizations of the WMAP best-fit cosmological model [26], only accounting for pure 
scalar perturbations {Cf^ = 0). We then introduce an anisotropic rotation of the 
Stokes parameters according to Eq. ( |T7j) . with the rotation angle given by Eq. (l39l) 
where we fix amax = 1-98° while {6'„, </>„} = {45°, 90°}. Angles are given in the Galactic 
coordinate system with 6 G [—90°, +90°] and G [0°,360°]. Furthermore, note that the 
value chosen for a^nax is consistent with presently available upper limits on the isotropic 
version of the birefringence model [H [261 123 US]- In figure [1] we report the map of the 
birefringence angle (see Eq. ( 1391) ) as input to the simulations. The map points out the 
symmetries expected for the proposed birefringence model, and in fact we have verified 
that, averaging over 1000 simulations, the full-sky power spectra of the rotated maps 
are consistent with those of the unrotated ones. 



Input a (e,0) = (45°,9O°) 




-2.0 ^_ 2.0 deg 



Figure 1. Map of the birefringence angle (see Eq. ([23)) as used in the simulations. 
The map has been generated in the Galactic coordinates, using the Healpix scheme 
[5S] with resolution parameter iVgido = 1024. 

We add to the rotated CMB maps described above isotropic white noise based on 
the Planck 143 GHz channel sensitivity [29], assuming a 30 months long mission. We 
choose this particular frequency channel because it is the one with the highest signal to 
noise ratio in polarization at high resolution. By following the procedure described in 
the previous section, we apply our set of masks to the simulated maps, computing for 
each mask the power spectra and deriving the best-fit of ao by minimizing the in E<1- 
f H3|) . In figure [2] the recovered values of are projected onto a map and, as one may 
see, they trace a smooth distribution with the highest values concentrated, as expected, 
around the input values of {9n,(t>n}- 

As described in the previous section, we can estimate the best-fit cimax from the 
maximum of this map. We derive the uncertainties on ctmax by means of 500 MC 
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Figure 2. Distribution on the sphere of the estimated ao. It has been obtained by 
analyzing CMB plus noise simulations with a birefringence angle ([M]) set by a^ax is 
1.98° and {6'„,0„} — {45°, 90°}. Note that here we are plotting the absolute value of 
the rotation angle ag ■ The map has been generated in the Galactic coordinates of the 
Healpix scheme [2S] with the resolution parameter Nside = 1024. 



simulations of unrotated CMB plus noise, the same simulations we use to assess the 
covariance matrices of the power spectra {Cf^). Note that considering unrotated CMB 
maps is just an approximation, nevertheless we expect that the contribution of the 
rotation itself to the power spectra covariances is subdominant. For the analyzed maps 
with the anisotropic rotation, we can individuate the disk i where the birefringence 
effect is maximum and compute the best-fit ao* for each of the MC simulations on that 
specific mask. Using these 500 values of aoi we can trace a "frequentist" probability 
distribution for the parameter and evaluate the one sigma error as the standard deviation 
of these values However, when we build the ao map as in figure El several, partially 
overlapping, disks contribute to constrain Oq, other than the particular disk i where the 
estimates is found to be maximal. Estimates from these disks contribute to lower the 
variance associated to ao with respect to erg, but they are obviously correlated due to 
overlap. In order to estimate the latter variance, we model the correlation degree in the 
overlapping disks as the fractional overlapping area between each pair of contributing 
disks. In other words, this correlation is if two disks are totally disjointed, and 1 in 
the limit of perfectly overlapping disks, all intermediate values being allowed. We then 
take as one sigma error: 

a = crs\/(W^C-iW)-i, (44) 

where W=[l,...,l] is a design matrix, while is the correlation matrix between 
OQi of different overlapping masks computed as described above. (Note how the rescaling 
factor reduces to the well familiar in the unrealistic case of N disjointed disks.) 

In figure [3] we present the histogram that traces the probability distribution of otmaxi as 
derived from the 500 MC simulations, and superimposed to it there is also the x^ia) 
curve obtained from the standard analysis of Eq. (H3ll on the considered maximum effect 



14 



mask (both rescaled by the factor a/ (W^C~^W)~i given in Eq. f l44p ) . The width of 
the histogram is compatible with that of the x^i^^) curve, confirming that neglecting 
the rotation in the MC simulations has a negligible impact. We want to stress that 
from figure [3] is evident that we can detect the birefringence effect (vertical red line) at 
high significance. Once the uncertainty on amax has been assessed, we can derive the 
confidence intervals on {9n, 4>n} (see figure Hj) as explained at the end of the previous 
section, i.e. by slicing the function ao{6, (p) in correspondence of the 1 and 2a errors on 
ao (which are related to the uncertainties on amax through Eq. (l40l) ). 

(e,0) = (45°,9O°) 




«max (deg) 

Figure 3. Histogram of ctmax estimated for 500 MC simulations of umotatcd CMB 
plus noise on the mask where the rotation effect is maximum. The solid curve in the 
plots represents the x^(q^) derived from the standard analysis on the same maximum 
effects mask, while the a in the label is the standard deviation of the values in the 
histogram. The plotted quantities have been rescaled by the factor ^ (W^C^^W)"! 
given in Eq. (|44p . The vertical line is the (imax evaluated from maps where the 
birefringence effect is present. 

Our main results are summarized in Table [H showing that the analysis recovers 
quite well the input parameters of the simulations. This means that our method of 
analysis is able to detect an anisotropic rotation effect (with amax given by the current 
upper limits on the isotropic case) and give reliable constraints on its parameters for 
data which are compatible with the Planck sensitivity. 

Oimax i ^CUm.ax 

1.97 ±0.20 [23.9,63.0] [49.5,114.7] 

Table 1. Estimates obtained for the three parameters of the model (all expressed in 
degrees). A0„ and A0„ have been derived at 2 cr of amax- We recall that the input 
value for Umax is 1.98° and {6'„,(/.„} = {45°, 90°}. 

In order to check for the robustness of the analysis, we perform further consistency 
tests. For instance, we verify that results remain stable by (reasonably) changing the 
range of multipoles taken into account in Eq. ( H3l) . In particular, results reported in 
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Table [Hand in the previous figures have been obtained by considering £ G [100,2000], 
where imin is bound by the fact that larger scales are not properly constrained by 
power spectra computed on small regions of the sky, while Imax corresponds roughly to 
the smallest angular scale that can be recovered at the 20% of the instrumental beam 
transfer function (which, for the Planck 143 GHz channel, we assume to be a Gaussian 
with FWHM=7.1' [29]). Moreover we derive constraints of the three parameters of the 
model by analyzing separately TB and EB rotated spectra. As foreseable, constraints 
obtained from EB spectra are tighter, since the cosmic variance of these spectra, which 
we recall have been computed on small patches of the sky, is smaller than that of TB 
spectra. Still for checking purposes, we have analysed the two hemispheres of the rotated 
maps separately, obtaining compatible estimates of the parameters. Results of all these 
checks are reported in Table El 





^max i ^^max 




A0„ 


TB 


2.38 ±0.34 


[21.9, 72.2] 


[33.1,137.8] 


EB 


1.83 ±0.22 


[25.2,70.7] 


[49.8,113.6] 


N 


2.25 ±0.28 


[22.3,66.3] 


[59.3,112.3] 


S 


1.89 ±0.29 


[25.2, 73.6] 


[27.5, 135.2] 



Table 2. Consistency checks on the estimated model parameters. The first two rows 
refer to the results obtained analysing Dj^ and Df^ spectra separately. The last two 
rows, instead, show results derived on a single hemisphere (respectively N=North and 
S=South). A0„ and A0„ have been derived at 2 ct of otmax- All values reported in 
this table are expressed in degrees. We recall that the input value for Umax is 1.98° 
and {On,4'n} = {45°, 90°}. The results of the full analysis are reported in Table 1. 
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A natural continuation of this analysis is to apply it to presently available data, 
such as those collected by WMAP, and we will address this point in a future paper. 

5. Conclusions 

It was recently shown that data from CMB polarization observations can be used 
to constrain quantum-gravity effects producing anomalous light propagation, with a 
sensitivity that is sufficient to test effects originating at the Planck scale. However, these 
observations were thought to be able to provide limits that are not really competitive 
with the ones coming from astrophysics, due mainly to the lower energy characterizing 
CMB photons. 

We established here that CMB data can also provide powerful constraints on non- 
isotropic effects, since they provide information from radiation coming from almost 
all directions on the sky. And in this case they can be truly competitive with the 
corresponding constraints on non-isotropic anomalous light propagation obtainable from 
observations of astrophysical sources, since in astrophysics one is limited to gaining 
information on only a few directions of propagation of signals (see also the discussion 
in [18]). 

For the specific non-isotropic anomalies for light propagation we considered here, 
which are expected in the quantum gravity- inspired models of Ref. [7l[18], we found that 
they don't produce any signature in the CMB full-sky power spectra, due to cancellations 
induced by the peculiar symmetries of the effect. Nevertheless, the anomalous effects can 
be anyway detected through a method that still exploits the power spectra formalism 
while operating on small patches of the sky in order avoid the cancellations. 

The effect we studied produces a rotation of the CMB polarization direction 
(birefringence), whose amount depends on the observation direction. We characterized 
the effect in terms of the amount of rotation in the direction of maximum effect (amax) 
and the two angles identifying this direction {9n and (pn)- These parameters are related 
in a simple way to the parameters characterizing the Lagrangian describing our model 
(see Eq. (l39l)). 

We have shown that data coming from the Planck satellite will be able to constrain 
anomalous non-isotropic light propagation with a sensitivity of 0.2 degrees on the 
amount of rotation in the direction of maximum effect and will also allow to identify 
the special direction pointed by the symmetry breaking vector with an uncertainty of 
roughly 40 degrees on the 6 angle and 60 degrees on the (j) angle. 

The sensitivity on a.max translates into a sensitivity on the parameter that indicates 
how far the constraint is from the Planck scale, i. e. the module of the symmetry-breaking 
vector n. We expect this parameter to be of order one if the effect is generated at the 
Planck scale. 

Exploiting the relation in Eq. ( l39i) . and taking into account also the correction due 
to photon redshift as described in [3], the sensitivity on Omax can be translated into a 
sensitivity on \n\ of the order of 10~^, which is one order of magnitude better than the 
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level required to test the Planck scale. 

Note also that the sensitivity can be further improved exploiting the availability of 
different energy channels in the Planck observations and the predicted dependence of 
the effect on the square of the photons energy. 
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